THE NUMBER OF LATIN SQUARES OF ORDER 11 



ALEXANDER HULPKE, PETTERI KASKI, AND PATRIC R. J. OSTERGARD 

Abstract. Constructive and nonconstructive techniques are employed to enu- 
merate Latin squares and related objects. It is established that there are 

(i) 2036029552582883134196099 main classes of Latin squares of order 11; 

(ii) 6108088657705958932053657 isomorphism classes of one-factorizations of 
^11,11; (iii) 12216177315369229261482540 isotopy classes of Latin squares of 
order 11; (iv) 1478157455158044452849321016 isomorphism classes of loops of 
order 11; and (v) 19464657391668924966791023043937578299025 isomorphism 
classes of quasigroups of order 11. The enumeration is constructive for the 
1151666641 main classes with an autoparatopy group of order at least 3. 



1. Introduction 

1.1. Latin squares. A Latin square of order n is an n x n array over n symbols 
such that every row and column is a permutation of the symbols. Latin squares 
are among the most fundamental combinatorial objects, and have many equivalent 
formulations jlj Theorem 1.11]. 

Latin squares are easily seen to exist for every n. It is, however, considerably 
more challenging to count the number of distinct Latin squares. In fact, the exact 
number is known only for n < 11, with the most recent case n = 11 settled by 
McKay and Wanless [H] . 

Here our interest is to enumerate Latin squares up to equivalence, whereby we 
focus on main classes (see j ]1.3[ ). We refer to |TS] for a brief historical account of 
enumerating small {n < 6) Latin squares; this history goes as far back as the 18th 
century and includes pioneering combinatorialists such as Cayley and Euler. There 
is a unique main class for n < 3, two main classes for n = 4, 5, and 12 main classes 
for n = 6. There are 147 main classes for n — 7; this classification was carried out by 
Norton [T5] and corrected by Sade [T7]. Kolesova, Lam, and Thiel [31 showed that 
there are 283657 main classes for n = 8. All these results were constructive, that 
is, representative Latin squares were (or could have been) obtained from each main 
class. Recently, McKay, Meynert, and Myrvold [13], on the other hand, counted 
in a nonconstructive way the number of main classes for n = 9 and n = 10; these 
numbers are 19270853541 and 34817397894749939, respectively. 

In this paper, we enumerate the main classes for n = 11 together with related 
objects such as quasigroups and loops (see [13] and 

Theorem 1. There are 

(i) 2036029552582883134196099 main classes of Latin squares of order 11; 

(ii) 6108088657705958932053657 isomorphism classes of one-factorizations of 
^11,11; 
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(iii) 12216177315369229261482540 isotopy classes of Latin squares of order 11; 

(iv) 1478157455158044452849321016 isomorphism classes of loops of order 11; 
and 

(v) 19464657391668924966791023043937578299025 isomorphism classes ofqua- 
sigroups of order 11. 

1.2. Group actions. We assume familiarity with group actions [11]. The Orbit- 
Stabilizer Theorem provides the methodological foundation for counting in |13j as 
well as in this paper; a similar approach is also used in [7 in the context of one- 
factorizations of a complete graph. 

Let r be a finite group that acts on a finite set . For an X e denote the 
orbit of X by [X] = {jX : 7 S T} and the stabihzer of X in F by Fx = {7 € T : 
'yX — X}. For a subgroup 11 < F, denote the conjugacy class of 11 in F by [U] = 
{7117^1 : 7 e r} and the normahzer of 11 in F by Fq = {7 G F : 7X17"^ = 11} . The 
Orbit-Stabilizer Theorem connects the size of an orbit with the size of a stabilizer 
subgroup by |F| = |[X]|.|Fx|. 

In the context of combinatorial enumeration, the elements of ^ are called "la- 
beled" objects, and the orbits of F on ^ are called "unlabeled" objects. Our 
interest is to count the unlabeled objects. To this end, denote by Ni the number 
of orbits [X] with |Fx| = i- Because ^ partitions into orbits, we have 

(1) lifi^in^^. 

i 

If we now know |F|, |^|, and Ni for each i > 2, then we can solve ([!]) for A^i 
and obtain the number of unlabeled objects J^i-^i- earlier studies (71 113j. 

the approach is to determine \ by a nonconstructive counting technique and the 
values Ni for i > 2 by a constructive classification up to isomorphism. 

In this paper, we extend the nonconstructive tools to the setting of prescribed 
symmetry, which enables us to determine the values Ni nonconstructively also for 
i > 1. There are two high-level ideas. First, the stabilizer subgroups of the elements 
of an orbit form a conjugacy class of subgroups of F. Thus, we can split the task 
of determining Ni into subtasks indexed by conjugacy classes [11] with |n| = i. 
Second, we can solve each such subtask via an analog of ([I]). Namely, instead of 
counting individual labeled objects X e ^ as in ([!]), we count formal pairs {X, E) 
consisting of a labeled object X E ^ and a subgroup S < Fx with E e [11]. 

Let us now proceed with a more detailed treatment. Denote by iVp] the number 
of orbits [X] with Tx G [H]. Taking the sum over all conjugacy classes [11] with 
[n[ = i, we have 

(2) ^» = E ^[n]- 

[n]:|n|=i 

For a conjugacy class [11], define the following two sets of formal pairs: 
if[n] = : X G if, S < Fx, S e [H]}, 

(3) 

%] = {{X, E) : X G if, S < Fx, S G [H]} . 

For X G if , let 



(4) 5[n]W = {S<Fx :Sg [n]}. 
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Observe that Sin]{jX) = 7S'[n](X)7"^ for all 7 e T. By the Orbit-Stabilizer 
Theorem, we thus have 

(5) \^[u]\-m2^ I J. I , F[n]|-|r| 2^ ir^i • 

[x] ' ^' [x]:|n|<|rx| ' ' 

Because Tx e [B] implies both S[u]iX) = {r^} and |rx| = |B|, we have (of. ([T 

(6) |^[n]h|r|^ + |^[n]|. 

Equation ^ is now the crux of our counting approach: First, we compute 
|.5f[n]| nonconstructively. Second, we constructively enumerate the orbits [X] with 
|B| < iTxl, which gives us |'^[n]| via Q and Finally, we solve Q for N^u]- 

One further simplification is that the task of computing \^[n] \ can be reduced 
from [B] to a representative group B < F in the conjugacy class. To this end, let 

(7) ^u = {X ■.U< Fx} . 

From (|3|, the fact that 7.ifn = ^yu-y-^ foi' all 7 e F, and the Orbit-Stabilizer 
Theorem, we have 

(8) |i?[n]l= E \^^\ = \m-\^u\^P-A^u\. 

We now proceed to apply this approach in the setting of Latin squares. 

1.3. Main classes. We work with the standard triple system representation of 
Latin squares. To this end, let R — {ri, r2, . . . , r„}, C = {ci, C2, . . . , c„}, and 
S = {si, S2, . . . , s„} be three pairwise disjoint n-element sets. The intuition is that 
R indexes the rows, C indexes the columns, and S is the set of symbols. We use 
design-theoretic terminology and call the elements of i? U C U S* points and the sets 
i?, C, S point classes. 

A Latin square of order n can now be represented as a set of triples, where 
each triple {ri,Cj,Sk} indicates that the symbol Sk occurs in row rj, column Cj of 
the array. Put otherwise, a Latin square (in triple system representation) is a set C 
of triples over RUCU S such that (a) any two points from different point classes 
occur together in a unique triple; and (b) no triple contains two points from the 
same point class. 

Assuming R, C, S are fixed but arbitrary, denote by ^ the set of all Latin squares 
in triple system representation. A "Latin square" in what follows refers to an 
element of Accordingly it is convenient to assume in what follows that a "pair" 
(of points) refers to a set of two points from two different point classes, and a 
"triple" (of points) refers to a set of three points from three different point classes. 

Let F be the group consisting of all permutations oi RU C U S that fix the 
partition {R, C, S} setwise. The structure of this group is a wreath product F = 
Sn I S3 = (Sn)'^ » S3 with the three copies of Sn permuting the points within 
each of the three point classes, and S3 permuting the point classes. In particular, 
|F| =3!(n!)3. 

The orbits of F on ^ are called main classes. Two Latin squares in the same 
main class are paratopic. The stabilizer of a Latin square in F is called the au- 
toparatopy group; the elements of the group are autoparatopisms. 
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1.4. Proof outline. The rest of this paper constitutes a proof of Theorem [T] Our 
approach is to determine iV^ for i — 2 via nonconstructive techniques enabled by j ]1.2| 
and for i > 2 via constructive enumeration of main class representatives. Finally, 
we solve for A^i using Q, with |Jf| given by [14j . 

We begin in Sj2] by establishing the necessary background on autoparatopisms 
of Latin squares. In f|3j we carry out a constructive enumeration of all main 
classes of Latin squares admitting an autoparatopy group of order at least 3. In 
^ we carry out a nonconstructive enumeration of Latin squares admitting fixed 
autoparatopisms of order 2. Finally, in fj5]we complete the proof of Theorem [l] by 
applying Q, Q, ([5|, and ([6]) to the results of ^}3]and ^ 

2. Groups of autoparatopisms 

2.1. Prime order. A group element a € F of prime order p is characterized up 
to conjugation in F by (a) the prime order p; (b) the number of points fixed in 
each of the point classes R, C, and S; and (c) the number of point classes R, C, S 
fixed. Denote by //j, fc, and fs the number of points fixed by a in R, C, and S, 
respectively. Up to conjugacy we may assume Jr > fc > fs- Denote by F the 
number of point classes R, C, S fixed by a. 

We proceed to narrow down the possible types {p, fn, fc, fs, F)- The results in 
the following two lemmata are well known; see, for example, [131 Theorem 1]. 

Lemma 2. The order p divides 3 — F. For each parameter f in fji, fc, fs it holds 
that / = or p divides n — f; furthermore, the latter property must hold for at least 
F of the parameters. 

Proof. A cycle decomposition of a consists only of p-cycles and fixed points. So if 
a point class is fixed — there are F such point classes — and contains / fixed points, 
then p divides n — f. Similarly, p divides S ~ F. □ 

Lemma 3. Any group of autoparatopisms of a Latin square that fixes points from 
at least two point classes has the property that the fixed points induce a subsquare. 
In particular, if fp > 1 and /c > 1 for a nonidentity autoparatopism, then //? = 
fc^fs< n/2. 

Proof. Consider the triples in a Latin square. If a triple has two fixed points relative 
to a group of autoparatopisms, then also the third point must be fixed; otherwise 
the two fixed points would occur in at least two triples, which is impossible. Thus, 
the number of triples with three fixed points is fnfc = fnfs = fcfs, from which 
we conclude fn — fc — fs when fR,fc > 1- Moreover, the triples with three fixed 
points induce a Latin square on the fixed points. The order of a proper subsquare 
in a Latin square of order n is at most n/2 [H Theorem 1.42]. □ 

We must consider the prime orders p = 2,3,5,7,11 for n — 11. Combining 
Lemmata [2] and |3j we can exclude all other types {p, fn, fc, fs,F) except those 
listed in Table [l] Table [T] also displays data related to the constructive enumeration 
to be discussed in fj3) 

2.2. Order four. We next determine the subgroups 11 < F of order 4 whose ele- 
ments of order 2 belong to the conjugacy classes in Table [T] We call such elements 
of order 2 eligible elements. 
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Type 


P 


/fl 


/c 


h 




T 


Seeds 


Mt 


1 


2 


1 


1 


1 


3 


- 


- 


- 


2 


2 


3 


3 


3 


3 


- 


- 


- 


3 


2 


5 


5 


5 


3 


- 


- 


- 


4 


2 


1 








1 


- 




- 


5 


rt 
I 


6 


(J 


(J 


1 








6 


2 


5 








1 








7 


2 


7 








1 








8 


2 


9 








1 








9 


2 


11 








1 








10 


3 


2 


2 


2 


3 


13 


13 


4 


11 


3 


5 


5 


5 


3 


13 


30 


25 


12 


3 














c2n 


339638 


165 


13 


5 


1 


1 


1 


3 


13 


7 


1 


14 


7 


4 


4 


4 


3 


13 


6 


16 


15 


11 











3 


ml 


28 


33 


16 


11 


11 








3 


flml 


1 


242 



We first compute the conjugacy classes of elements of F (there are 34048 classes) 
using standard algorithms and isolate the 374 classes of order 4. Amongst these we 
find 6 classes for which the squares of elements are eligible. We observe that the 
elements in each of these classes are conjugate to their inverses. Therefore there 
are 6 classes of cyclic groups of order 4. They are listed as types 38-43 in Table [3j 
with generator permutations (column "Generators" ) appearing in Table m 

We note that the conjugacy classes of eligible elements, as well as of elements of 
order 4 with eligible squares, are determined by their permutation cycle structure 
l°i2°2 4<»4 for oi + 2a2 + ia^ = 3n and the number F of point classes R, C, S that 
they fix. Table [2] indicates these 16 conjugacy classes. 

To classify elementary abelian subgroups of order 4, we note that each such 
subgroup must contain three elements of order 2, any two of which generate the 
group. We will first classify such generating pairs (a,/?), where both a G T and 
/3 S r are eligible, and (3 is contained in the centralizer of a, which we denote by 
Fc, = {7 e r : 707"^ — a}. Furthermore a/3 must be eligible as well. Clearly, 
any two of the three elements a, /3, a/3 generate the same group. We therefore can 
assume without loss of generality that (in an arbitrary total order on the conjugacy 
classes of F, for example, the "Type" number in Table [2]) the class of a is not larger 
than the class of /3, which is not larger than the class of a/3. If this is fulfilled, we 
call the pair (a, /3) valid. 

Furthermore, we need to classify such pairs only up to conjugacy in F. For 
this, we may assume that a is the chosen representative of its conjugacy class and 
/3 is chosen up to conjugacy by F^. We determine the representatives for /3 by 
computing for each eligible a the conjugacy classes of Fq, and identifying amongst 
these the classes of eligible elements. For each such class representative /3 we test 
whether (a,/3) is valid, and collect the valid pairs in a list. 
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Table 2. Autoparatopism types for order dividing 4 



Type 


Cycles 


F 


1 


j^32l5 


3 


2 


1^212 


3 


3 


J^1529 


3 


4 


^l2l6 


1 


r 

5 


1 3ol5 
12 


1 


6 


15214 


1 


7 


1^213 


1 


8 


1^21^ 


1 


9 


111211 


1 


10 


2^32346 


3 


11 


1^2246 


3 


12 


1^2146 


3 


13 


^946 


3 


14 


II2446 


1 


15 


2^32346 


1 


16 


^33 


3 



Because conjugacy of subgroups conjugates elements to elements, the lexico- 
graphic choice of generating pairs ensures that the only remaining potential conju- 
gacy amongst subgroups generated by such pairs can happen if both a and /3, or /3 
and aP are in the same class under F. 

A classification of such pairs on the computer finds 37 valid pairs; none of the 
corresponding subgroups are conjugate, as verified by explicit conjugacy tests in 
r. These groups are listed as types 1-37 in Table [3j again generators appear in 
Table m 

Calculations were done in the system GAP [3]. 

Up to conjugacy in F, there are therefore exactly 43 groups H < F of order 4 
whose elements of order 2 belong to the conjugacy classes in Table [T] 

It turns out that these groups are distinguished up to conjugacy in F by the 
conjugacy classes of their elements (column "Elements") and the lengths \°-i2°-^4°"^ 
for ai -I- 2a2 -I- 4a4 = 3?! of their orbits on i? U C U 5' (column "Orbits"). 

We can rule out the following groups with combinatorial arguments. 

Lemma 4. The groups of types 1, 3, 4, 5, 9, 10, 13, 16, 24, 27, 40, and 41 are 
not admitted by a Latin square of order 11. 

Proof Let i? = {1, 2, . . . , 11}, C = {12, 13, ... , 22}, and S = {23, 24, ... , 33}. 

Types 4, 10, 16, 27, 40, and 41: All point classes contain fixed points, but their 
number varies between point classes. This contradicts Lemma [3j 

Type 1: Let 

i?i = {1, 2, 3, 4, 5, 6} , i?2 = {7, 8, 9, 10, 11} , 

Ci = {12, 13, 14, 15, 16, 17} , C2 = {18, 19, 20, 21, 22} , 

Si = {23, 24, 25, 26, 27, 28} , ^2 = {29, 30, 31, 32, 33} . 

Lemma [3] applied to 71 reveals that i?2 U C2 U 6*2 induces a 5 x 5 subsquare. Thus, 
there exists no triple with exactly two points from i?2 U C2 U 6*2. 
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Table 3. Candidate groups of order 4 



7 



Type 


Generators 


Elements 


Orbits 


T 


Seeds 


Mt 


1* 


71.72 


3, 6, 9, 


16 


l'"2"4'^ 


- 


- 


- 


2 


71.73 


3, 4, 7, 


16 


-l^lglO^S 


flm2 


102 


11 


3* 


71.74 


2, 3, 3, 


16 


l'^2 


- 


- 


- 


4* 


71.75 


2, 3, 3, 


16 




- 


- 


- 


5* 


71.76 


3, 7, 8, 


16 




- 


- 


- 


6 


71.77 


3, 5, 6, 


16 




flm2 


230 


11 


7 


71.78 


1, 2, 3, 


16 




13 


40 


1 


8 


71.79 


3, 6, 7, 


16 


1^2'' 4'' 


f3n 


80 


10 


9* 


71 . 710 


3, 4, 5, 


16 


■1^12844 


- 


- 


- 


10* 


71.711 


2, 3, 3, 


16 


-1^521241 


- 


- 


- 


11 


7i . 7i2 


3, 5, 8, 


16 


l^2^A^ 


f3n 


80 


1 


12 


7i . 7i3 


3, 6, 7, 


16 


l'^2 4 


f3n 


80 


1 


13* 


71.714 


2, 3, 3, 


16 


1^2^4^ 


- 


- 


- 


14 


71.715 


3, 5, 6, 


16 


r2^4:^ 


f3n 


80 


1 


15 


7l6. 717 


2, 2, 2, 


16 


1^4'' 


13 


12 


9 


16* 


7l6. 718 


2, 2, 2, 


16 


■1^52644 


- 


- 


- 


17 


7l6. 7l9 


1. 1. 2, 


16 


l'^2'^4° 


13 


20 


1 


18 


7l6. 720 


2, 5, 9, 


16 


-^32744 


f3n 


127 


1 


19 


716,721 


2, 7, 7, 


16 


1^2'' 4'' 


f3n 


83 


1 


20 


7l6. 722 


2, 5, 5, 


16 


1^2^4'^ 


f3n 


83 


1 


21 


7l6. 723 


2, 6, 8, 


16 


^32744 


f3n 


127 


1 


22 


716.724 


2, 4, 8, 


16 


^^12844 


flm2 


608 


11 


23 


7l6. 725 


2, 6, 6, 


16 


1^2*4* 


flml 


2215 


22 


24* 


7l6, 726 


2, 4, 4, 


16 


1^2*4^ 


- 


- 


- 


25 


7l6. 727 


2, 5, 7, 


16 


1^2*4'' 


flml 


2835 


22 


26 


7l6. 728 


2, 2, 2, 


16 


1^2® 4^ 


13 


10 


1 


27* 


7l6. 729 


2, 2, 2, 


16 


l''2^4^ 


- 


- 


- 


28 


716.730 


2, 5, 7, 


16 


132^4^ 


f3n 


127 


1 


29 


7l6. 731 


2, 6, 6, 


16 


1^2^4^ 


f3n 


83 


1 


30 


7l6. 732 


2, 4, 6, 


16 


1^2^4^ 


flml 


2854 


22 


31 


716.733 


2, 5, 5, 


16 


12 4 


flml 


3351 


22 


32 


734 . 735 


1, 4, 9, 


16 


■,^12645 


flml 


194 


22 


33 


734, 736 


1, 6, 7, 


16 


■,^12645 


flml 


2300 


22 


34 


734.737 


1, 5, 8, 


16 


^^12645 


flml 


1083 


22 


35 


734 , 738 


1. 4, 7, 


16 


1^2*4« 


flml 


1696 


22 


36 


734, 739 


1, 5, 6, 


16 


^12446 


flml 


6406 


22 


37 


734,740 


1, 4, 5, 


16 


■^12247 


flml 


1761 


22 


38 


741 


2, 13, 13, 


16 


^946 


13 


20 


9 


39 


742 


2, 10, 10, 


16 


■1^32346 


f3 


20 


1 


40* 


743 


2, 11, 11, 


16 


1^224^ 








41* 


744 


2, 12, 12, 


16 


^72146 








42 


745 


2, 15, 15, 


16 


-1^32346 


f3n 


4 


1 


43 


746 


2, 14, 14, 


16 


^12446 


flml 


242 


22 



*) Excluded by Lemma [i] 
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Table 4. Generators for Table |3] 



R = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11} 

C = {12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22} 

S = {23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33} 



71 = (1,2)(3,4)(5,B){12,13)(14,15)(16,17){23,24)(25,26)(27,28) 

72 = (12,23)(13,24){14,25)(15,26)(16,27){17,28)(18,31)(19,29){20,30)(21,32)(22,33) 

73 = (8,10)(9,11)(12,23){13,24)(14,25)(15,26)(16,27)(17,28)(18,31){19,29)(20,30)(21,32)(22,33) 

74 = (1,2)(8,10)(9,11)(16,17)(19,22)(20,21)(27,28)(29,33)(30,32) 

75 = (1,2)(8,10)(9,11)(12,14)(13,15)(19,22)(23,25)(24,26)(29,33) 

76 = (1,2)(12,23)(13,24){14,25)(15,26)(16,27)(17,28)(18,31)(19,29){20,30)(21,32)(22,33) 

77 = (1,2)(8,10)(9,11)(12,23)(13,24)(14,25)(15,26)(16,27)(17,28)(18,31)(19,29){20,30)(21,32)(22,33) 

78 = (3,5)(4,6)(8,10)(9,11)(12,14)(13,15){19,22)(20,21)(23,25){24,26)(29,33)(30,32) 

79 = (3,5)(4,6)(12,23)(13,24)(14,25)(15,26)(16,27)(17,28)(18,31)(19,29)(20,30){21,32)(22,33) 

710 = (1,2)(3,5)(4,6){8,10)(9,11)(12,23)(13,24)(14,25)(15,26)(16,27){17,28)(18,31)(19,29)(20,30)(21,32)(22,33) 

711 = (1,2)(8,10)(9,11)(16,17)(19,22)(20,21)(23,25)(24,26)(29,33) 

712 = (9,11)(12,23){13,24)(14,25)(15,26)(16,27)(17,28)(18,31)(19,29){20,30)(21,32)(22,33) 

713 = (1,2)(9,11)(12,23)(13,24){14,25)(15,26)(16,27){17,28)(18,31)(19,29){20,30)(21,32)(22,33) 

714 = (3,5)(4,6)(9,11)(12,14)(13,15){19,22)(23,25)(24,26){29,33) 

715 = (3,5)(4,6)(9,11)(12,23)(13,24){14,25)(15,26)(16,27)(17,28)(18,31)(19,29){20,30)(21,32)(22,33) 

716 = (1,2)(3,4)(5,6){7,S){12,13)(14,15)(16,17)(18,19)(23,24)(25,26){27,28)(29,30) 

717 = (1,3)(2,4)(5,7){6,8){12,14)(13,15)(16,18){17,19)(23,25)(24,26){27,29)(28,30) 

718 = (1,3)(2,4)(5,7){6,8){12,15)(13,14)(16,17){21,22)(23,26)(24,25){27,28)(32,33) 

719 = (1,3)(2,4)(5,7){6,8){10,11)(12,14)(13,15){16,18)(17,19)(21,22)(23,25)(24,26)(27,29)(28,30)(32,33) 

720 = (12,23)(13,24)(14,25)(15,26)(16,27){17,28)(18,29)(19,30)(20,31)(21,32)(22,33) 

721 = (5,6)(7,8)(12,23)(13,24)(14,25)(15,26)(16,27)(17,28)(18,29)(19,30)(20,31){21,32)(22,33) 

722 = (1,3)(2,4)(5,7){6,8){12,23)(13,24)(14,25){15,26)(16,27)(17,28){18,29)(19,30)(20,31)(21,32)(22,33) 

723 = (7,8)(12,24)(13,23){14,26)(15,25)(16,28){17,27)(18,29)(19,30){20,31)(21,32)(22,33) 

724 = (10,11)(12,23){13,24)(14,25)(15,26)(16,27)(17,28)(18,29){19,30)(20,31)(21,33)(22,32) 

725 = (5,6)(7,8)(10,11)(12,23)(13,24)(14,25)(15,26)(16,27)(17,28)(18,29)(19,30){20,31)(21,33)(22,32) 

726 = (1,3)(2,4)(5,7){6,8){10,11)(12,23)(13,24){14,25)(15,26)(16,27){17,28)(18,29)(19,30)(20,31)(21,33)(22,32) 

727 = (7,8)(10,11)(12,24){13,23)(14,26)(15,25){16,28)(17,27)(18,29){19,30)(20,31)(21,33)(22,32) 

728 = (1,3)(2,4)(7,8){10,11)(12,15)(13,14)(16,17)(21,22)(23,25){24,26)(29,30)(32,33) 

729 = (1,3)(2,4)(7,8){10,11)(12,14)(13,15)(16,18)(17,19)(23,26)(24,25)(27,30)(28,29) 

730 = (1,3)(2,4)(12,25)(13,26)(14,23)(15,24)(16,27)(17,28)(18,29)(19,30)(20,31){21,32)(22,33) 

731 = (1,3)(2,4)(7,8){12,26)(13,25)(14,24){15,23)(16,28)(17,27){18,29)(19,30)(20,31){21,32)(22,33) 

732 = (1,3)(2,4)(10,11)(12,25)(13,26)(14,23)(15,24)(16,27)(17,28)(18,29)(19,30){20,31)(21,33)(22,32) 

733 = (1,3)(2,4)(7,8){10,11)(12,26)(13,25)(14,24)(15,23)(16,28)(17,27)(18,29)(19,30){20,31)(21,33)(22,32) 

734 = (1,2)(3,4)(5,6){7,8){9,10){12,13)(14,15)(16,17){18,19)(20,21)(23,24)(25,26)(27,28)(29,30){31,32) 

735 = (12,23)(13,24)(14,25)(15,26)(16,29){17,30)(18,27)(19,28){20,31)(21,32)(22,33) 

736 = (7,8)(9,10)(12,23)(13,24)(14,25)(15,26)(16,29)(17,30)(18,27)(19,28){20,31)(21,32)(22,33) 

737 = (3,4)(12,23)(13,24)(14,25)(15,26)(16,29){17,30)(18,27)(19,28){20,31)(21,32)(22,33) 

738 = (7,9)(8,10)(12,23)(13,24){14,25)(15,26)(16,29){17,30)(18,27)(19,28){20,31)(21,32)(22,33) 

739 = (3,4)(7,9)(8,10)(12,23)(13,24){14,25)(15,26)(16,29){17,30)(18,27)(19,28){20,31)(21,32)(22,33) 

740 = (1,2)(3,5)(4,6){7,9){8,10){12,23)(13,24)(14,25){15,26)(16,29)(17,30){18,27)(19,28)(20,31){21,32)(22,33) 

741 = (1,2,3,4)(5,6,7,8)(12,13,14,15){16,17,18,19)(23,24,25,26)(27,28,29,30) 

742 = (1,2,3,4)(5,6,7,8){9,10){12,13,14,15)(16,17,18,19)(20,21)(23,24,25,26)(27,28,29,30)(31,32) 

743 = (1,2,3,4)(5,6,7,8){12,13,14,15){16,17,18,19)(20,21)(23,24,25,26)(27,28,29,30)(31,32) 

744 = (1,2,3,4)(5,6,7,8){12,13,14,15){16,17,18,19)(23,24,25,26)(27,28,29,30)(31,32) 

745 = (1,13,2,12)(3,15,4,14)(5,17,6,16)(7,19,8,18)(9,20)(10,21)(11,22)(23,24,25,26){27,28,29,30) 

746 = (1,13,2,12)(3,15,4,14)(5,17,6,16)(7,19,8,18)(9,20)(10,21)(11,22){23,24,25,26){27,28,29,30){31,32) 



LATIN SQUARES OF ORDER 11 



9 



Next consider pairs with one point from Ci and one from Si ; there are |Ci | • |5i | — 
36 such pairs. The action of the group partitions these pairs combined with any 
points from R into orbits, exactly six of which have length 4; for example, one such 
orbit is 

{{x, 12, 25}, {x, 14, 23}, {71(2;), 13, 26}, {71(0;), 15, 24}} . 

As 71 (x) ~ X when x € R2, any such x can be in at most one orbit (otherwise it 
would occur in at least 2-4 = 8 triples together with points from Ci but |Ci| = 
6 < 8). 

It follows that x £ R2 in at most five of the above mentioned six orbits of length 
4, whereby x e i?i in at least one of the orbits. Such an orbit with x G Ri contains 
two pairs {x, y} with y € Ci (and recall that the third point of the triple is from 
^i). Finally, consider how pairs {x, z} with z E S2 are covered. The third point in 
the corresponding triple has to come from Ci (by the initial comment). However, 
we have 5 points in S2 to be combined with x and 6 — 2 = 4 points in Ci , which is 
not possible. 

Type 3: It follows by Lemma[3]that there are 5x5 subsquares over {7, 8, 9, 10, 11}, 
{18, 19, 20, 21, 22}, {29, 30, 31, 32, 33} and {3, 4, 5, 6, 7}, {12, 13, 14, 15, 18}, {23, 24, 
25,26,31} for 71 and 74, respectively. 

Consider the four triples containing the point 3 and one point from {19,20, 
21, 22}. By the given subsquare structure, the only possibilities for the third point 
of the triple are the points {27,28}, implying that {3,27} or {3,28} would occur 
in more than one triple. 

Type 5: The proof for Type 1 is applicable as such, since the only difference 
between these cases is the additional transposition (1,2) in generator 75, which 
leads to different orbits^for example, {1, 12, 25}, {2, 14, 23}, {2, 13, 26}, {1, 15, 24} 
instead of {1, 12, 25}, {1, 14, 23}, {2, 13, 26}, {2, 15, 24}— but has no impact on the 
arguments of the proof. 

Type 9: The only point fixed by 710 is 7, which is then the only possible third 
point in triples containing the points in transpositions of 710 with one point each 
from C and S. Since there are 11 such triples, the point 7 does not occur in any 
other triples. 

By considering possible choices for x in {x, 12, 24}, we get up to symmetry four 
cases with the orbits 

{{1,12, 24}, {2, 13, 23}}, 

{{3, 12, 24}, {4, 13, 23}, {5, 13, 23}, {6, 12, 24}} , 

{{8, 12, 24}, {8, 13, 23}, {10, 13, 23}, {10, 12, 24}} , and 

{{7, 12, 24}, {7, 13, 23}}. 

Because of the earlier argument, we cannot have x = 7, and because pairs of points 
occur twice in two other cases, we must have x G {1, 2}. 

Repeated use of the previous reasoning gives that each of the pairs {12,24}, 
{13,23}, {14,26}, {15,25}, {16,28}, and {17,27} must occur in triples together 
with the points 1 and 2 (three times each). Now consider triples {l,x,y} with 
y £ {29,30,31,32,33} = ^2. It follows from the subsquare implied by 71 that we 
cannot have a; e {18, 19, 20, 21, 22} = C2. There are therefore |C \ C2I = 6 possible 
choices for x, 3 of which already occur in triples together with the point 1. But 
6-3 = 3<5 = 1521. 
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Type 13: By applying Lemma [3] to 71 and 71714 we get that there is a 5 x 5 
subsquare over {7,8,9,10,11}, {18,19,20,21,22}, {29,30,31,32,33} and a 3 x 3 
subsquare over {7,8, 10}, {18,20,21}, {30,31,32}, respectively. But the latter is a 
subsquare of the former, and a 5 x 5 square cannot contain a 3 x 3 subsquare [1] 
Theorem 1.42]. 

Type 24: Analogously to the case of Type 9, we argue that, as the only point 
fixed by 726 is 9, this point occurs exactly in triples with the other two points from 
transpositions of 725. 

By considering possible choices for x in {x, 12, 24}, we get up to symmetry three 
cases with the orbits 

{{1, 12, 24}, {3, 13, 23}, {4, 12, 24}, {2, 13, 23}} , 

{{10, 12, 24}, {10, 13, 23}, {11, 13, 23}, {11, 12, 24}} , and 

{{9, 12, 24}, {9, 13, 23}}. 

But because of the earlier argument, we cannot have x = 9, and pairs of points 
occur twice in the other two cases. □ 

We remark that the combinatorial proofs in Lemma[4]were triggered by computer 
searches showing nonexistence. 

3. Constructive enumeration 

3.1. Admissible groups. We recall the following corollary of the Sylow Theorems. 

Theorem 5. |16t Corollary 4.15] Let K he a finite group, let p he a prime, and let 
k he a positive integer. If p^ divides \K\, then A contains a suhgroup of order p^ . 

Let us say that a Latin square C admits A < L as a group of autoparatopisms if 
A < r£ . Combining Theorem [5] with the analysis in ^|2j we have that every Latin 
square C of order 11 with \Tc\ > 3 admits at least one subgroup conjugate to (a) 
a group of odd prime order in Table [T] or (b) a group of order 4 in Table |3] not 
excluded by Lemma |4j Let us call these groups the admissihle groups. 

3.2. The enumeration. We apply the framework in [5 to construct exactly one 
representative from each main class of Latin squares admitting at least one admis- 
sible group. 

Associated with each conjugacy class [11] of admissible groups is a set of seeds, 
at least one of which is guaranteed to occur in every Latin square admitting 11. 

The formal definition of a seed is as follows. Let H < F be an admissible group. 
Let T C i? U C U 5 be a set of points whose size and composition is determined by 
column "T" in Tables [T] and |3] as follows: 

ii : indicates i points fixed by 11; 

Tui : indicates i points moved by 11; 

ci : indicates i points in same point class; 

n : see Item (d) below. 

(For example, "flm2" indicates that T consists of 3 points, 1 point fixed by H, and 
2 points moved by Ft.) Finally, let 5 be a union of H-orbits of triples such that, 
referring to the elements of S as hlocks, 

(a) any pair occurs in at most one block; 

(b) each point in T occurs in exactly n blocks; 
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(c) T has nonempty intersection with at least one block on every H-orbit on 
S] and 

(d) the set T occurs in at least one block unless the composition of T has the 
"n" indicator. 

Each such tuple (11, T, 5) is a seed associated with the conjugacy class [11]. 

Let 7 e r act on a seed by 7(n, T, S) = (707^^, 7T, 7^). The orbits of T on the 
set of all seeds are the isomorphism classes of seeds. A Latin square C contains (or 
extends) a seed (H, T, 5) if H < and S C £. 

Given a seed {Il,T,S), the task of finding all Latin squares that extend the 
seed is an instance of the exact cover problem. In particular, the task is to cover 
exactly once all uncovered pairs using H-orbits of triples; each triple covers the 
pairs occurring in it. 

Our constructive enumeration approach now proceeds as follows. First, we clas- 
sify the seeds up to isomorphism using the algorithms described in [5] . The number 
of nonisomorphic seeds associated with each conjugacy class of admissible groups 
is given in column "Seeds" in Tables [l] and |3] Once the seeds have been classified, 
we use libexact |8j to find all extensions of each seed to a Latin square. As each 
extension £ of a seed (H, T, S) is constructed, it is subjected to isomorph rejection 
tests developed in [S], which derive from the canonical augmentation technique of 
McKay [T^. First, we identify a canonical F£-orbit of seeds contained in C. If 
(n, r, S) does not occur in the identified orbit, we reject C from further considera- 
tion. Second, we test whether C is the lexicographic minimum of its F(n.T.5)-orbit. 
If not, we reject C from further consideration. Otherwise we accept C as the unique 
representative of its main class. See [5] for a detailed exposition of classification of 
combinatorial objects. 

From an implementation perspective the present algorithm is almost identical 
to the one used in [3 [7j ■ Some additional implementation effort was required to 
work with the groups of order 4, and to implement associated correctness checks 
(to be discussed in j |3.3[ ). In particular, we must list all subgroups of order 4 in 
Tc- This was implemented essentially using brute force, that is, by iterating over 
all elements of order 4, and all pairs of elements of order 2 that generate a group 



of order 4. (This strategy suffices because Fc is in most cases small; see Table 12 ) 
In the search we find 105670178597 Latin squares that extend the seeds. Among 

these, we find 1151666641 main classes. The number of main classes Ni for each 

i > 3 is given later in Table [12] 

The search was distributed to a network of 155 Linux PCs using the batch system 

autoson [TT]. In total the search consumed about 1.3 years of CPU time. 



3.3. Correctness. We carry out a double counting check to gain confidence in the 
correctness of the constructive enumeration. The validation procedures considered 
in this section and later in |4.4| are essential as, for example, number theoretic 
properties discussed in (TB] and elsewhere are not applicable here. 

Fix a conjugacy class of admissible groups [11]. We compute in two different 
ways the total number of tuples ((E, T, 5), £) where (E,T, 5) is a seed contained 
in C and satisfying Eg [11]. 

For a Latin square C, denote by seedsjn] (>C) the number of seeds associated with 
[n] contained in £. For a seed (E,T, 5), denote by ext(E,T, 5) the number of 
extensions of the seed into a Latin square. By the Orbit-Stabilizer Theorem, we 
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Table 5. Double count values for odd prime orders 



Type 


P 


Jr 


fc 


fs 


F 


Count 


10 


3 


2 


2 


2 


3 


88699187523260511795806208000000 


11 


3 


5 


5 


5 


3 


12309174115893617098752000000000 


12 


3 














19601984696323546934786654208000000 


13 


5 


1 


1 


1 


3 


710224896233056606617600000 


14 


7 


4 


4 


4 


3 


1186258276475008450560000000 


15 


11 











3 


5968708870624483737600000 


16 


11 


11 








3 


38160882055721779200000 



have 

[C] ' ^' [(S,T,5)]:Se[n] ' 

where the sum on the left-hand side is over the main classes of Latin squares, and 
the sum on the right-hand side is over the isomorphism classes of seeds associated 

with [n]. 

The right-hand side of ^ is accumulated for each classified seed (n,T, 5). In 
particular, ext(n,T, 5) is simply the number of solutions found in the exact cover 
search when extending (H, T, S). 

The left-hand side of ^ is accumulated whenever a constructed C is accepted 
as the representative of its main class; that is, for every [IT] we compute seeds[n]('C) 
and accumulate accordingly. To compute seeds [n] (C), we iterate over the subgroups 
S < Fc with S S [n]. For each such E, we accumulate seeds[n]('C) by the value 
listed in column M-r in Tables [T] and |3] To justify this, first observe that whenever 
S and T are fixed, C determines 5 in a seed (E, T, S) occurring in C. In particular, 
S is the union of all S-orbits that contain a triple that has nonempty intersection 
with T. Now, whenever S is fixed, Mt counts the number of eligible sets T. The 
values Mt can be determined using combinatorial arguments based on the structure 
of n and the size and composition constraints for T in relation to 11. 

Let us give two examples to illustrate the combinatorial arguments. First, Type 
12 in Table [I] has Mt — 165 because there are 3(2^) = 165 ways to select two 
distinct points in a common point class. Second, Type 15 in Table |3] has Mt — 9 
because T must consist of 3 points fixed by E and must occur as a subset of (that 
is, must be equal to) a triple. There are 9 choices for such a T by Lemma[3j namely 
the 9 triples of the subsquare of order 3 induced by the points fixed by E. 

The computed left-hand and right-hand sides of ([9| agree for each [11] in our 
constructive enumeration, which gives us confidence that the results are correct. 
We display the double count values for reference in Tables [5] and [6] 

4. NONCONSTRUCTIVE ENUMERATION 

4.1. One-factorizations with symmetry. Our objective in this section is to 
reduce the task of computing \^u\ for subgroups 11 < F with |n| = 2 to the task 
of counting one-factorizations of a complete bipartite graph with forced symmetry. 

We refer to [19] for basic graph-theoretic terminology. All graphs considered are 
undirected, loopless, and without parallel edges. A one-factor of a graph G is a 
spanning 1-regular subgraph of G. A one-factorization of G is a set of one-factors 
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Type Count 



2 


671631524180703313920000000 


6 


2014894572542109941760000000 


7 


3516906890255319171072000000 


8 


54951670160239362048000000000 


11 


1831722338674645401600000000 


12 


5495167016023936204800000000 


14 


21980668064095744819200000000 


15 


^ f\f^ r* r\ A Ar\r'f^^r\'\r\r\Airit^A a a r\r\r\r\r\r\ 

307619449557019948744704000000 


17 


39732092772896030588928000000 


18 


13737917540059840512000000000 


19 


268164150381968086794240000000 


20 


65282584150364362113024000000 


21 


179508789190115249356800000000 


22 


16790788104517582848000000000 


23 


1663496959090765967917056000000 


25 


501843074867821516161024000000 


26 


18539743730111373901824000000 


ZQ 


OOODoOo4o ( zoo4zd44Uo2DooUUUUUU 


29 


568713151711703904288768000000 


30 


246354443069481975545856000000 


31 


573841974259992911413248000000 


32 


23802621216964125445324800000 


33 


22878993544910822248022016000000 


34 


1477320700587875009298432000000 


35 


464500362123374411907072000000 


36 


12596852214923599074557952000000 


37 


61521447614952423555072000000 


38 


919890958482406920683520000000 


39 


39728022278810086932480000000 


42 


45133638424943262695424000000 


43 


2532319498770923774803968000000 



of G such that every edge of G occurs in a unique one-factor in J^. Denote by LF{G) 
the number of distinct one-factorizations of G. 

Let Kn,n be the complete bipartite graph with vertex set RU C and bipartition 
{R,G}. Let be a one-factorization of Kn.n, and let / : S* — > be a bijection 
that "labels" the one-factors in T with elements of S. 

Lemma 6. There is a one-to-one correspondence between the tuples , /) and the 
Latin squares in . 

Proof. Each triple {r^, Cj, Sfe} in a Latin square L corresponds to the edge {r^, Cj} 
in the one-factor f{sk) in T . □ 

In particular, Lemma |6] implies 

1^1 = n! • LF(i^„,„) . 

We now introduce forced symmetry into this setting, that is, we proceed to study 
in light of Lemma |6] Our task is fortuitously simplified by the assumption 
|n| = 2, which implies that 11 fixes at least one of the point classes R, C, S setwise. 
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Without loss of generality (up to conjugation of 11 in F), in what follows we will 
assume that 11 fixes S setwise. 

Let $ be the group consisting of all permutations oi RUC that fix the partition 
{R, C} setwise. Let $ act on the set of spanning subgraphs of Kn.n by permuting 
the vertices, and extend the action in the natural way to one-factorizations of such 
graphs. 

Observe that each element tt e H can be restricted to i? U C because 11 fixes S 
setwise. Let A be the restriction of H to i? U C, and observe that A < <i>. 

Consider now a Latin square C with 11 < Fc, and let (J^, /) be the corresponding 
tuple given by Lemma [6] We clearly have A < $jr. However, it is not the case 
that any tuple (J^, /) with A < <i>jr corresponds to a Latin square with H < r£. In 
particular, the action of A on need not be "compatible" with the action of 11 on 
S. We proceed with a detailed analysis. 

Let be a one-factorization of Kn^n with A < <i>jF, and let / : S* ^ J-" be a 
bijection. For any S G A, define 6 as the permutation oi RU C U J-' induced by the 
action oi S on RU C and on J^. (We assume that RU C and T are disjoint sets.) 
Note that 5 is well-defined because S e A < Let A = {(5 : 5 G A}. Extend 
/ by the identity mapping on i? U C to the mapping / : RU C U S ^ RUC U F. 
Let f-^Af = {/"iJ/ -.S e A}, and observe that /-^A/ < F. We say that (J", /) 
agrees with FT if f^^Af — H. 

Lemma 7. The tuple {J-, /) agrees with FE if and only if the corresponding Latin 
square C satisfies 11 < F^ . 

Proof. Follows from Lemma [6] because n fixes S' setwise. □ 

We say that T is compatible with 11 if there exists an f : S ^ T such that (J-", /) 
agrees with FI. Denote by ^ the group of all permutations oi RU C U S that fix 
RU C pointwise. 

Lemma 8. Let J- be compatible with II. Then there are exactly |^'n| bijections f 
.such that (J-, f) agrees with 11. 

Proof. Fix a reference bijection and establish a one-to-one correspondence between 
the set of all bijections and "^n- D 

Denote by LF(_ftr„^„, FT) the number of distinct one-factorizations of Kn.n com- 
patible with n. From Lemmata [7] and [8] it thus follows that 

(10) l^nl = |*n| •LF(if„,„,n). 

The next three lemmata rely on our assumption that 11 has order 2. 

Lemma 9. Suppose H restricted to S has oi orbits of length 1 and 02 orbits of 
length 2. Then 

(11) |*n| =ai!a2!2°^. 

Proof. There are a\\a2\'2!^'^ permutations in '5 that fix the nonidentity element of 
n under conjugation. □ 

Lemma 10. Let T be a one- factorization of Kn.n with A < ^jr. Then J- is 
compatible with H if and only if the action of A on J-^ has an equal number of 
orbits of each length as the action of H on S . 
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Proof. The permutation of S induced by the nonidentity element of 11 can be rela- 
belled (by conjugation with an appropriate /) to the permutation of induced by 
the nonidentity element of A if and only if there are an equal number of cycles of 
each length. □ 

Let G be a spanning subgraph of Kn^n- We say that a one- factorization F of G 
is a (A, a)-factorization if (a) A < $jr; and (b) for each i = 1,2, the action of A on 
J-" has Ui orbits of length i with a — {ai, 02). Denote by LF(G, A, a) the number of 
distinct (A, a)-factorizations of G. 

Lemma 11. Suppose II restricted to S has ai orbits of length 1 and ai orbits of 
length 2. Then 

(12) LF(i^„,„,n) = LF(i^„,„,A,a). 

Proof. Immediate from Lemma [TO} □ 



We have now reduced the task of computing |^n| via ( [10[ ), (11), and (12 1 to 
computing LF(if„,„, A, a). 

We employ two different methods to compute the values LF(if„_„, A, a) for the 
conjugacy classes of order 2 in Table [T] The first method stems from a recursion 
over regular spanning subgraphs of Kn^n that has evolved over a period of more 
than half a century into the form presented by McKay and Wanless [13]. The 
second method extends a "forward accumulation" technique used in [7j to count 
the number of distinct one-factorizations of K14. 

Our present contribution is to modify both methods to take into account the 
symmetry forced by A and the required structure a = (01,02) for the lengths of 
the A-orbits of one-factors. 



4.2. Backward recursion. For convenience in what follows, let us abbreviate 
6*1 = (1,0) and 6*2 = (0,1). We say that a spanning subgraph G of -ftr„ „ is 
(A, a)-factorizable if G has at least one (A, a)-factorization. Observe that a (A, a)- 
factorizable graph is necessarily fc-regular with k — '^^ia.i. Furthermore, each 
(A, a)-factorizable graph decomposes into an edge-disjoint union consisting of oi 
(A, ei )-factorizable and 02 (A, e2)-factorizable graphs. Denote by I{G, A, i) the set 
of all (A, ei)-factorizable spanning subgraphs of G. 

Counting in two different ways the number of distinct (A, a)-factorizations of G 
with one individualized A-orbit consisting of i one-factors, we obtain, for > 1, 

(13) •LF(G,A,a) = ^ LF(F, A, e^) • LF(G - F, A, a - e,) . 

FeI(G,A,i) 

In particular, if we know (for example, by recursion) the value LF(_ff, A, a — e.i) 



for all (A,a — ei)-factorizable H, we can use (13) to compute LF(G, A,a) for all 
(A, a)-factorizable G. Furthermore, LF(G,A,a) needs to be computed for only 
one representative (G, A) chosen from the $-orbit [(G, A)]. Here $ acts on G by 
permutation of vertices and on A by conjugation. 

This is analogous to the method (without the forced A and a) used in [T3] and 
other earlier studies on counting Latin squares. 
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4.3. Forward accumulation. We can visit every (A, a)-factorizable orbit [(G, A)] 
with cii > 1 by the following procedure: for each (A,a — ei)-factorizable orbit 
[(iJ, A)], consider exactly one tuple {H,A) from the orbit; for each such tuple 
{H, A), consider each graph F e I{Kn.n — H, A, i); for each tuple {H, F, A), visit 
the (A, a)-factorizable orbit [{H U F, A)]. 

To compute the value LF(G, A, a) for each visited orbit [(G, A)], we associate 
with [(G, A)] an accumulator variable a;[((3 a)] that is initially set to zero and in- 
cremented whenever [(G, A)] is visited. Our objective is to have the value • 
LF(G, A, a) in X[(g,A)] when the procedure terminates. 

We proceed to analyze the visiting procedure, with the objective of determin- 
ing an appropriate increment to the counter variable on each visit. The follow- 
ing two lemmata are immediate consequences of the Orbit-Stabilizer Theorem (cf. 
[21 Lemmata 1,2]). 

Lemma 12. Any tuple (G, A) in the orbit [{H U F, A)] admits exactly 

a(i7,F,A)-'*(^"^'^)' 



\^{H,F,A)\ 

decompositions G = H' U F' into tuples {H' , F' , A) in the orbit [(iJ, F, A)]. 
Lemma 13. The procedure visits an orbit [{H U F^ A)] exactly 

T{H,F,A)- 



times via tuples (H' , F' , A) in the orbit [{H, F, A)]. 



It now follows from Lemma 12 and ( |13[ ) that, for a.; > 1, 
(14) 

a, • LF(G, A, a) = ^ LF(i^, A, e,) • LF(i/, A, a - e,) 

(H,F,A):HUF=G 

J2 a{H,F,A)-LF{F,A,e,)-LF{H,A,a~e,). 

[{H,F,A)]:HUF=G 

This observation enables us to accumulate the value Oi ■ LF(G, A, a) to x\^(^g,a)]- 
Namely, each time [(G, A)] is visited via a tuple {H, F, A), we increment x^cA)] 
the rule 

(15) x^^GA)] ^ H(GA)] + 1[h'fa) ■ ^' "^'^ ■ ^' " "^'^ ■ 

Equivalently, for each tuple (H, F, A) considered by the visiting procedure, we apply 
the rule 

\^{HUFA)\ 



(16) X[(^HuF.A)] ^ X[(^HuF,A)] + "HS '~T~ ' ^^(F, A, Ci) • LF(i7, A, a - Ci) . 

l*(-fi",A)l 

Analogously to [71 Lemma 3], we have the following result. 
Lemma 14. The total accumulation to a)] «s Oi ■ LF(G, A, a). 



Proof. By Lemma 13 and (15), the total accumulation to X[(^ga)] fro™ orbit 



[{H,F,A)] satisfying G ^ H U F is a{H,F,A) ■ LF{F,A,e^) ■ LF{H, A,a - Ci). 



Taking the sum over all such classes, the claim follows by (14). □ 
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4.4. Implementation and correctness. Associated with each graph F in §4.2| 



and S4.3 there are exactly LF(i^, A, e^) distinct A-orbits (of length i) of one-factors; 
that is, the distinct (A, ei)-factorizations of F. From an implementation perspective 
it is convenient to iterate over these A-orbits directly instead of the underlying 



graphs F. In particular, we achieve the required accumulation in (131 and ( 16 ) 
without computing the values LF(_F, A, e;) explicitly. 

Starting with the empty spanning subgraph of K^^m we use the forward accumu- 
lation method to compute the values LF(G', A, a) for all (A, a)-factorizable orbits 
[(G, A)]. We carry out the forward accumulation in two different ways. First, we 
accumulate the A-orbits in order of increasing length; that is, we first add all A- 
orbits of length 1, and then the A-orbits of length 2. Second, as a consistency 
check, we add the orbits in order of decreasing length. 

In addition to the forward accumulation method, we verify the results using 



the recursion (13), that is, when we have accumulated the value LF(G', A,a) by 



addition of orbits of length i, we apply (13 1 to check the result using the values 
LF(ff,A,a-e*). 

The two different orders of adding orbits enable yet another way to cross-check 
the results (cf. [5] and [71). Namely, 

- A)] ui^i ■ LF(G, A, a) ■ LF(X„,„ - G, A, b) 
(17) LF(X„^„,A,a + fe)= 

holds for every nonnegative a = (01,02) and b = (61,62) with '^^i{ai + bi) = n, 
where the sum is over all (A, a)-factorizable orbits [(G, A)]. Due to the two different 
orders of adding orbits, we cannot utilise all possible decompositions a + b for 
checking, however. For example, if we add orbits in the orders 1,1,1,2,2,2,2 and 



2,2,2,2,1,1,1, then we can use eight different decompositions in (17). 

The approaches were implemented using the C programming language with the 
help of three software libraries: nauty [IQi for isomorphism testing, the GNU Mul- 
tiple Precision Arithmetic Library [4 for handling large integers (and intermediate 
rationals), and libexact [& for simplifying the search for A-orbits of one-factors. 
The number of orbits [(G, A)] at any step of the search never exceeded 200000, so 
memory requirement was not a major issue of the search (compare with [71 113)). 
The computations took about 20 days on a contemporary Linux PC. 

4.5. Results. The numbers of (A, a)-factorizable orbits [(G, A)] are shown in Ta- 
bles [7] and jS] for situations when the one-factor orbits are appended in order of 
increasing and decreasing length i, respectively. The leftmost column contains the 
number of one-factor orbits a^, and there is one column for each of the conjugacy 
classes of order 2 in Table [TJ 

The second column of Tabic [lOj summarizes the results of the search, which 
enable us to determine the values |-Sf[n] I in the rightmost column of Table 
the combination of (jsj), (jlOj), (jlll), and (IT2I. That is, 



10 



via 



\^m\ - 7|5^-ai!a2!2'^^ •LF(ifii,ii,A,(ai,a2)), 
|ln| 

where |r| = 3!(11!)^ and the values jPnl are displayed in Tablejoj these are derived 
by considering the number of permutations in P that fix the nonidcntity element of 
n under conjugation. 
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Table 7. Number of (A, a)-factorizable orbits, increasing length 





\ lype 1 


o 
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6 
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10 
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22 


22 


22 


22 
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16 
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1365 
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30429 


30429 


30429 
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20 
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46 
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146085 


146085 
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11 
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49 
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125763 


125763 


7 
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45 
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19567 


19567 


8 






1 
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36 
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630 
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21 


21 


10 
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1 


11 


















1 



Table 8. Number of (A, a)-factorizable orbits, decreasing length 



Ei^AType 1 


2 


3 


4 


5 


6 


7 


8 


9 
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1 
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1 
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1 
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6 
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6 
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6 


6 
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1464 


3 


266 


266 


266 


266 


1089 


22 


3 


34792 


1182 


1 
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170321 
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1365 


4 
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16 


5 


6 


6 


33851 


187719 
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5 
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10 


10 


1 


46 
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6 


1 


1 
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1 


5 
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Table 9. |rn| for conjugacy classes with |n| = 2 



Type 






irni 


1 


(1!- 5! -25)3. 


3! 


339738624000 


2 


(3! • 4! • 2^)3 . 


3! 


73383542784 


3 


(5! • 3! • 23)3 . 


3! 


1146617856000 


4 


1!- 5! •25. 11! 


• 2 


306561024000 


5 


3! • 4! • 2* • 11! 


• 2 


183936614400 


6 


5! -3! -23 -11! 


• 2 


459841536000 


7 


7! • 2! • 22 • 11! 


• 2 


= 3218890752000 


8 


9! • 1! • 2^ • 11! 


• 2 


= 57940033536000 


9 


ll!-0!-2°- 11! 


• 2 


= 3186701844480000 
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Table 10. |^[n]| for conjugacy classes with |n| = 2 
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Type 



LF(Xn,ii,A,a) 



l-^[n]l 



1 157811617463135109120 

2 16482860057738870784 

3 57076088832000000 

4 81867724734136320 

5 91614680894074060800 

6 984580848455978188800 

7 588014650532826316800 

8 39448015149028147200 

9 252282619805368320 



680681075601465561779007809126400000 
197485300889612684060962848768000000 
109414881030165485322240000000000 
391332044120479692845639270400000 
437923008879888637579625496576000000 
4706348408721523581497414516736000000 
2810741057151330357801067216896000000 
188563593954724468142515421184000000 
1205924234796020705924638310400000 



Table 11. I'^pjl and A^p] for conjugacy classes with |n| — 2 



Type 


l^[n]| 




1 


1788612130139530711950950400000 


3567419044431 


2 


2090765465832065147338752000000 


1035003382971 


3 


40157841013697699905536000000 


573229534 


4 


37761778205327560959590400000 


2050761205 


5 


1253232559784675856089088000000 


2295134347303 


6 


3238510153755989660663808000000 


24665809156818 


7 


2028400471919270993854464000000 


14730996126592 


8 


252164628166875895037952000000 


988254629517 


9 


14823428114332518344294400000 


6320133846 


N2 




47291560812217 



5. Number of equivalence classes 

5.1. Main classes. As a result of the computations described in ^ and J|4j we 
have a constructive enumeration for all main classes with an autoparatopy group of 



order at least 3. The numbers of such main classes are given in Table 12 for i > 3 



Table 10 enables us to solve the number of main classes with an autoparatopy group 
of order 2. Indeed, we compute |^[n]| by direct summation over the constructed 
main classes ([5]), and use ^ together with Table 10 to arrive at Table 
Table Inland Table [l2] (for i > 3) give 



11 



(18) 



iriE 



= 9023583561995938862980803959193600000 . 



i>2 



(19) 



The number of distinct Latin squares of order 11 is [14] 

1^1 = 11! • 10! • 5363937773277371298119673540771840 

= 776966836171770144107444346734230682311065600000 . 

Solving for A^i in ([T]) using (18) and (19), we obtain 

Ni = 2036029552535590421717241 
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Table 12. Main classes of Latin squares of order 11 



i 


Ni 


i 


Ni 


1 


2036029552535590421717241 


24 


274 


2 


47291560812217 


27 


27 


3 


1111651266 


30 
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39004721 


32 


16 
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9131 


36 


86 
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960771 
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1294 
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18 
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63 
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12 


4191 


72 


14 


14 
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108 


5 


15 


4 


110 


2 


16 


631 


120 


1 


18 


625 


162 


1 


20 


37 


324 


1 


21 


37 


7260 


1 


22 


4 


E.N. 


2036029552582883134196099 



and hence 

^ iV, = 2036029552582883134196099 . 

i 

Table [T2| gives the number Ni of main classes for each order i of the autoparatopy 
group. 

5.2. One-factorizations and isotopy classes. For a Latin square £, denote by 
Par(£) the autoparatopy group of C, and denote by Is(£) the subgroup of Par(£) 
that fixes each of the point classes R, C, S setwise. This group is called the au- 
totopy group of C, and its elements are called autotopisms. Let Ty(£) = 3, 2, 1, 1 
when |Par(£)|/|Is(£)| = 1,2,3,6, respectively. Observe that these quantities are 
independent of the main class representative C because Is(£) is a normal subgroup 
of Par(/:). 

By [13 Theorem 2(ii) and 2(iii)], each main class [£] splits into Ty(£) isomor- 



phism classes of one-factorizations of Kn,n (cf. S4.1) and into 3! • |Is(£)|/|Par(£)| 
isotopy classes. 

Taking the sum over all main classes, we obtain Theoremjljii) andjljiii). Table[l3] 
gives the order of the autoparatopy group and the order of the autotopy group for 
each main class. 

5.3. Quasigroups and loops. A permutation with cycles of length i for i > 1 
is said to have cycle structure {ai,a2, ■ ■ ■)■ For an autotopism a € Is(£), define 
/(q^) = rii^i'*"' ^^'^ 5(0^) = fli if Q; is not the identity permutation and has 
the same cycle structure (ai, 02, . . .) in each of the point classes R, C, S; otherwise 
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Table 13. Orders of autoparatopy and autotopy groups 
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The number of isomorphism classes of quasigroups is, by J13', Theorem 4], 

[C] ' ^ " aels(C) 

The number of isomorphism classes of loops is, by 13, Theorem 5], 



We obtain Theorem [l|iv) and[ljv) by (19), (20), (21 1, and summation over the 
main classes. 
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